Phase transitions of the generalized contact process with two absorbing states 
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We investigate the generalized contact process with two absorbing states in one space dimension 
by means of large-scale Monte-Carlo simulations. Treating the creation rate of active sites between 
inactive domains as an independent parameter leads to a rich phase diagram. In addition to the 
conventional active and inactive phases we find a parameter region where the simple contact process 
is inactive, but an infinitesimal creation rate at the boundary between inactive domains is sufficient 
to take the system into the active phase. Thus, the generalized contact process has two different 
phase transition lines. The point separating them shares some characteristics with a multicritical 
point. We also study in detail the critical behaviors of these transitions and their universality. 
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I. INTRODUCTION 

Many systems in physics, chemistry, and biology are 
far from thermal equilibrium, even if they are in time- 
independent steady states. In recent years, continu- 
ous phase transitions between different nonequilibrium 
steady states have attracted lots of attention. Just as in 
equilibrium, these transitions are characterized by large- 
scale fluctuations and collective behavior over large dis- 
tances and long times. Examples can be found, e.g., in 
surface growth, granular flow, chemical reactions, popu- 
lation dynamics, and even in traffic jams [H-0I- 

Continuous nonequilibrium phase transitions can be 
divided into different universality classes according to 
their critical behavior, and considerable effort has been 
devoted to categorizing the variety of known transitions. 
A well-studied type of nonequilibrium phase transitions 
separates fluctuating (active) steady states from absorb- 
ing (inactive) states where fluctuations stop completely. 
The generic universality class for these so-called absorb- 
ing state transitions is directed percolation (DP) Q. 
More specifically^ it was conjectured by Janssen and 
Grassberger [H, |T(| that all absorbing state transitions 
with a scalar order parameter and short-range interac- 
tions belong to this class as long as there are no extra 
symmetries or conservation laws. While nonequilibrium 
transitions in the DP universality class are ubiquitous in 
both theory and computer simulations, experimental ver- 
ifications were only found rather recently in ferrofluidic 
spikes fllj ] and in the transition between two turbulent 
states in a liquid crystal fl^ . 

Absorbing state transitions in universality classes dif- 
ferent from DP can occur in the presence of additional 
symmetries or conservation laws. Hinrichsen [l3| intro- 
duced nonequilibrium lattice models with n > 2 absorb- 
ing states. In the case of two symmetric absorbing states 
(n = 2), he found the transition to be in a new univer- 
sality class, the ^-symmetric directed percolation class 
(DP2). If the symmetry between the absorbing states is 
broken, the critical behavior reverts back to DP. In one 
dimension, the DP2 universality class coincides 0] with 
the parity-conserving PC class [14[ which is observed, 



e.g., in the branching-annihilating random walk with an 
even number of offspring (BARWE) [la ]. 

In this paper, we revisit one of the stochastic lattice 
models introduced in Ref. [ljj, the generalized contact 
process with two absorbing states in one space dimen- 
sion. Compared to the simple contact process [l6j], this 
model contains an additional dynamical process, viz., the 
creation of active sites at the boundary between domains 
of different inactive states. By treating the rate for this 
process as an independent parameter we uncover a rich 
phase diagram with two different types of phase transi- 
tions, separated by a special point that shares many char- 
acteristics with a multicritical point. We perform large- 
scale Monte-Carlo simulations of this model to study in 
detail the critical behavior of these transitions. 

Our paper is organized as follows. We introduce the 
generalized contact process with several absorbing states 
in Sec. [ill ln Sec. IIII1 we summarize the mean-field the- 
ory for this system. Sec. IIVI is devoted to the results 
and interpretation of our Monte-Carlo simulations. We 
conclude in Sec. W\ 



II. THE GENERALIZED CONTACT PROCESS 
WITH SEVERAL ABSORBING STATES 

The contact process [lj| is a paradigmatic model in the 
DP universality class. It is defined on a d-dimensional 
hypercubic lattice. Each lattice site r can be in one of 
two states, namely A, the active (infected) state or I, the 
inactive (healthy) state. Over the course of the time evo- 
lution, active sites can infect their nearest neighbors, or 
they can become inactive spontaneously. More precisely 
the contact process is a continuous-time Markov process 
during which active sites turn inactive at a rate /x, while 
inactive sites become infected at a rate Am/ (2c?) where 
m is the number of active nearest neighbors. The healing 
rate \x and the infection rate A are external parameters 
whose ratio determines the behavior of the system. 

If \i 3> A, healing dominates over infection. All in- 
fected sites will eventually become inactive, leaving the 
absorbing state without any active sites the only steady 
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state. Thus, the system is in the inactive phase. In the 
opposite limit, the infection survives for infinite 

times, i.e., there is a steady state with a nonzero density 
of active sites. This is the active phase. The nonequi- 
librium phase transition between these two phases at a 
critical value of the ratio \j \x is in the DP universality 
class. 

In 1997, Hinrichsen [l3[ introduced a generalization of 
the contact process. Each lattice site can now be in one 
of n+ 1 states, the active state A or one of the n different 
inactive states Ik (k = 1 . . . n). k is sometimes called the 
"color" index. The dynamics of the generalized contact 
process is defined via the following rates for transitions 
of pairs of nearest-neighbor sites, 

w(AA Alfe) = w(AA IfeA) = p,/n , (1) 

w(AI fe Ifelfe) = w(I fe A Ifelfe) = Hk , (2) 

w{Al k AA) = w{l k A -)• AA) = A , (3) 

w{l h li -> I fe A) = w(l k h -> Ah) = a , (4) 

with k, I = 1 . . . n and k ^ Z. All other rates vanish. We 
are mostly interested in the fully symmetric case, /ife = /i 
for all fc. For n = 1 and fl — fj,, the so defined generalized 
contact process coincides with the simple contact process 
discussed above. One of the rates p,, fi, A, and a can be set 
to unity without loss of generality, thereby fixing the unit 
of time. We choose A = 1 in the following. Moreover, to 
keep the parameter space manageable, we focus on the 
case p, = [i in the bulk of the paper. The changes for 
fi 7^ fi will be briefly discussed in Sec. |V] 

The process @ prevents inactive domains of different 
color (different k) to stick together indefinitely. They can 
separate, leaving active sites in between. Thus, this tran- 
sition allows the domain walls to move through space. It 
is important to realize that without the process ((4]), i.e., 
for a = 0, the color of the inactive sites becomes unim- 
portant, and all Ife can be identified. Consequently, for 
(7 = 0, the dynamics of the generalized contact process 
reduces to that of the simple contact process for all values 
of n. 

Hinrichsen [l3| studied the one-dimensional general- 
ized contact process by means of Monte-Carlo simula- 
tions, focusing on the case <r = A = 1. For n = 2, he 
found a nonequilibrium phase transition at a finite value 
of fi which separates the active and inactive phases. The 
critical behavior of this transition coincides with that of 
the PC universality class. For n > 3, he found the model 
to be always in the active phase. The Monte-Carlo simu- 
lations were later confirmed by means of a non-hermitian 
density-matrix renormalization group study [l7j . 

Motivated by a seeming discrepancy between these re- 
sults and simulations that we performed during our s tudy 
of absorbing state transitions on a percolating lattice [l8| , 
we revisit the one-dimensional generalized contact pro- 
cess with two inactive states. In contrast to the earlier 
works we treat the rate a of the process (j4{ as an inde- 
pendent parameter (rather than fixing it at a = A = 1). 



III. MEAN-FIELD THEORY 

To get a rough overview over the behavior of the gen- 
eralized contact process with two inactive states, we first 
perform a mean-field analysis. Denoting the probabilities 
for a site to be in state A, Ii, and I2 with P A , Pi, and 
Pi, respectively, the mean- field equations read: 

dP A /dt = (1 - n)P A -P 2 A + 2aP l P 2 , (5) 
dPx/dt = fiP A /2 - P A P X - aP x P 2 , (6) 
dP 2 /dt = nP A /2 - P A P 2 - oP x P 2 . (7) 

Let us begin by discussing the steady states which are 
given by the fixed points of the mean-field equations. 
There are two trivial, inactive fixed points P\ = 1, P A — 
P 2 = and P 2 = 1, P A = Pi = 0. They exist for all 
values of the parameters fi and a and correspond to the 
two absorbing states. In the case of a = 0, these fixed 
points are unstable for \i < 1 and stable for fi > 1. In 
contrast, for a > 0, they are always unstable. 

The active fixed point is given by P\ = P 2 and fulfills 
the equation 

= (1 - fi)P A - Pi + <r(l - P A ) 2 I2 . (8) 

For (7 = 0, this equation reduces to the well-known mean- 
field equation of the simple contact process, = (1 — 
/J,)P A — P A with the solution P A = 1 — /i for [i < 1. 
Thus, for a = 0, the nonequilibrium phase transition of 
the generalized contact process occurs at /i = tf? = 1. 
This means, it coincides with the transition of the simple 
contact process, in agreement with the general arguments 
given in Sec. [TTJ In the general case, a 7^ 0, the steady 
state density of active sites, P A , is given by the positive 
solution of 

P A = -— -- (l-fi-a± vV - 2m + 1 + 2m<t) . (9) 

We are particularly interested in the behavior of P A for 
small a. As long as fj, < /i£ p = 1 (i.e., in the active phase 
of the simple contact process), a small, nonzero a only 
provides a subleading correction to P A . At /i = [if — 1, 
the density of active sites vanishes as P A ~ ^/a with 
(7 — y 0. Finally, for fi > [if = 1, the density of active 
sites vanishes as P A ~ cr/(/i — 1). 

We thus conclude that within mean-field theory, the 
generalized contact process with two inactive states is 
in the active phase for any nonzero a. This agrees with 
older mean-field results but disagrees with more sophisti- 
cated methods which pred ict a nonequilibrium transition 
at a finite value of \i |l3l . \vK s . The mean- field dynamics 
can be worked out in a similar fashion. We find that 
the approach to the stationary state is exponential in 
time anywhere in parameter space except for the critical 
point of the simple contact process at [i — 1, a = 0. How- 
ever, it is known that mean-field theory does not reflect 
the correct long-time dynamics of t he g eneralized contact 
process which is of power-law type }13j . Therefore, we do 
not analyze the mean-field dynamics in detail. 
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IV. MONTE CARLO SIMULATIONS 
A. Method and overview 

We now turn to the main part of the paper, viz., large- 
scale Monte-Carlo simulations of the one-dimensional 
generalized contact process with two inactive states. We 
perform two different types of calculations: (i) decay runs 
and (ii) spreading runs. Decay runs start from a com- 
pletely active lattice; we monitor the time evolution of 
the density p(t) of active sites as well as the densities 
pi(t) and p2{t) of sites in inactive states Ii and I2, re- 
spectively Spreading simulations start from a single ac- 
tive (seed) site embedded in a system of sites in state Ii . 
(From a domain wall point of view, the spreading runs 
are therefore in the even parity sector.) Here we mea- 
sure the survival probability P s (t), the number of sites 
in the active cloud N s (t) and the mean-square radius of 
this cloud, R 2 (t). 

In each case, the simulation proceeds as a sequence of 
events. In each event, a pair of nearest-neighbor sites is 
randomly selected from the active region. For the spread- 
ing simulations, the active region initially consists of the 
seed site and its neighbors; it is updated in the course of 
the simulation according to the actual size of the active 
cluster. For the decay runs, the active region comprises 
the entire sample. The selected pair than undergoes one 
of the possible transitions according to eqs. ([1]) to dU 
with probability tw. Here the time step r is a constant 
which we have fixed at 1/2. The time increment associ- 
ated with the event is r /N pa i r where N pa i r is the number 
of nearest-neighbor pairs in the active region. 

Using this method we studied systems with sizes up 
to L = 10 6 lattice sites and times up to t max = 10 8 , 
exploring the parameter space < p < 1 and < a < 1. 
The a — p phase diagram resulting from our simulations 
is displayed in Fig. [TJ This phase diagram shows that 
the crossover from DP critical behavior at a = to DP2 
(or, equivalently, PC) critical behavior at a > occurs in 
an unusual fashion. The phase boundary cr c (p) between 
the active and inactive phases does not terminate at the 
critical point of the simple contact process located at 
(p,a) = (p c /,0) « (0.30325,0). Instead, it ends at the 
point (p, <j) = (p*,0) ~ (0.552,0). In the parameter 
range p cp < p < p* , the system is inactive at a = 0, but 
an infinitesimally small nonzero a takes it to the active 
phase. 

Thus, the one-dimensional generalized contact process 
with two inactive states has two types of phase transi- 
tions, (i) the generic transition occurring at p > p* and 
a = u c (p) > (marked by the dashed blue line and 
arrows in Fig. [T]) and (ii) the transition occurring for 
p cp < p < p* as a approaches zero (solid red line and 
arrows) . We note in passing that our critical healin g ra te 
for a = 1 is p c = 0.628(1), in agreement with Ref. [13[ 

In the following subsections we first discuss in detail 
the simulations that lead to this phase diagram, and then 
we present results on the critical behavior of both tran- 



1 

0.8 
0.6 
0.4 
0.2 





[■ 0.555 0.565 0.575 0.585 



Active 



Inactive 



0.1 



0.3 



0.5 



0.7 



FIG. 1: (Color online) Phase diagram of the ID general- 
ized contact process as function of the healing rate p and 
the boundary rate a. A line of DP2 (PC) transitions (blue 
dashed line) separates the active and inactive phases. For 
a — ► 0, this line does not terminate in the simple contact pro- 
cess critical point at p c J" ~ 0.30325 and but at p* « 0.552. 
For Pc P < p < p* , the system is inactive at a = (thick solid 
red line), but an infinitesimal a takes it to the active phase. 
Inset: Close to the endpoint at p* , the phase boundary be- 
haves roughly as a c ~ (p — p*) 2 . 



sitions as well as special point (p*,0) ~ (0.552,0) that 
separates them. 



B. Establishing the phase diagram 

We first performed a number of spreading simulations 
at d — and various p for maximum times up to 3 x 10 4 . 
The resulting number N s (t) of active sites in the clus- 
ter is shown in Fig. [5] The figure demonstrates that 
the transition between the active and inactive phases oc- 
curs at p = 0.30325(25). A fit of the critical curve to 
N s ~ i e =p yields 6 cp = 0.315(5). As expected from 
the general arguments in Sec. [Til both the critical heal- 
ing rate and the initial slip exponent Q cp agree very well 
with the results of the simple contact process (see, e.g., 
Ref. for accurate estimates of the DP exponents). 
Thus, at cr = 0, the generalized contact process under- 
goes a transition in the directed percolation universality 
class at p = pf = 0.30325(25). 

We now turn to nonzero a. Because the domain 
boundary process Q creates extra active sites, it is clear 
that the phase boundary between the active and inac- 
tive phases has to shift to larger healing rates p with in- 
creasing a. In the simplest crossover scenario, the phase 
boundary a c (p) would behave as a c ~ {p-p c c p ) 1 /'t' where 
is a crossover exponent. To test this scenario, we per- 
formed spreading simulations for times up to 10 7 at sev- 
eral fixed p > p cp in which we vary a to locate the tran- 
sition. Examples of the resulting N s (t) curves for several 
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FIG. 2: (Color online) Spreading simulations at a — 0: Num- 
ber 7V S of active sites as a function of time t. The solid 
line for /i = 0.30325 represents a fit to 7V S ~ t @cp yielding 
6 C p = 0.315(5). The data are averages over 25000 runs. 
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FIG. 3: (Color online) Spreading simulations: Number N 3 
of active sites as a function of time t for several a at fixed 
fi — 0.428 (panel a) and /i = 0.6 (panel b). The data are 
averages over 10 3 (at the smallest a) to 10 5 runs. 



a at /j, — 0.428 and fj, = 0.6 are shown in Fig. [3] The 
set of curves for /i = 0.6 (Fig. behaves as expected: 
Initially, N s (t) follows the behavior of the simple con- 
tact process at this /i. At later times, the curves with 
a 0.25 curve upwards implying that the system is in 
the active phase. The curves for a i; 0.25 curve down- 
ward, indicating that the system is in the inactive phase. 
Thus, u c {[i = 0.6) w 0.25. ' 
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FIG. 4: (Color online) Spreading simulations at a — 0.1 for 
several \i close to the phase boundary. Main panel: Number 
N s of active sites as a function of time t. Inset: Survival 
probability P s as a function of time t. The data are averages 
over 10 runs. 

In contrast, the set of curves for /i = 0.428 (Fig. [3ji) 
behaves very differently. After an initial decay, N s (t) 
curves strongly upwards for all values of a down to the 
smallest value studied, a — 10 -4 . This suggests that at 
fi = 0.428, any nonzero a takes the generalized contact 
process to the active phase. The phase transition thus 
occurs at a = 0. 

We determined analogous sets of curves for many dif- 
ferent values of the healing rate in the interval /i^. p = 
0.30325 < n < 0.65. We found that the phase transi- 
tion to the active phase occurs at a = for fi^ p < [i < 
fx* = 0.552, while it occurs at a nonzero a for healing 
rates fj, > /i*. This establishes the phase diagram shown 
in Fig. [1] The phase boundary thus does not follow the 
simple crossover scenario outlined above. In the follow- 
ing subsections, we analyze in detail the critical behavior 
of the different nonequilibrium phase transitions. 

C. Generic transition 

We first consider the generic transition occurring at 
fi > fj,* « 0.552 and nonzero a (the blue dashed line in 
Fig.[T]). Figure @] shows a set of spreading simulations at 
ct = 0.1 and several /j, in the vicinity of the phase bound- 
ary. The data indicate a critical point at /i m 0.582. We 
performed analogous simulations for several points on the 
phase boundary. Figure [S] shows the survival probability 
P s and number N s of active sites as functions of time for 
all the respective critical points. In log-log representa- 
tion, the N s and P s curves for different a and /i are per- 
fectly parallel, i.e., they represent power-laws with the 
same exponent. Fits of the asymptotic long-time behav- 
ior to P s = B a t~ s and N s = Gat® give estimates of 
5 = 0.289(5) and 8 = 0.000(5). Moreover, we measured 
(not shown) the mean-square radius R 2 (t) of the active 
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FIG. 6: (Color online) Critical density decay simulations: 
Density pa of active sites as function of t for several points 
(fj., a) on the generic phase boundary. The solid lines are fits 
to a power law pA = B„t~ a giving a — 0.285(5). The data 
represent averages of 400 runs with system size L = 10 4 . 



FIG. 5: (Color online) Critical spreading simulations: Sur- 
vival probability P a and number of active sites N 3 as func- 
tions of t for several points (/i, <r) located on the generic phase 
boundary. The inset shows the prefactor B a of the critical 
power law P a = B a t~ s as a function of a. The solid line is a 
fit to B a ~ <7~ c which gives £ = 0.284. 



versality class (see, e.g., Ref. 0, Q). We therefore con- 
clude that the critical behavior of the generic transition 
of generalized contact process with two inactive states is 
universally in this class. 



cloud as a function of time. Its long time behavior fol- 
lows a universal power law. Fitting to R 2 {t) ~ t 2 / z gives 
2/z = 1.145(5) (z = 1.747(7)). Here z = J /v x is the dy- 
namical exponent, i.e., the ratio between the correlation 
time exponent v\\ and the correlation length exponent v±. 

In addition to the spreading simulations, we also per- 
formed density decay simulations for several (/i, a) points 
on the phase boundary. Characteristic results are pre- 
sented in Fig. [Bl The figure shows that the density pa 
of active sites at criticality follows a universal power law, 
Pa = B a t~ a at long times. The corresponding fits give 
a = 0.285(5) which agrees (within the error bars) with 
our value of the survival probability exponent 6. We thus 
conclude that the generic transition of our system is char- 
acterized by three independent exponents (for instance 
v±, z and 5) rather than four (as could be expected for 
a general absorbing state transition [4j). We point out, 
however, that even though P s and pa show the same 
power-law time dependence at criticality, the behavior 
of the pref actors differs. Specifically, the prefactor B a 
of the density is increasing with increasing a while the 
prefactor B a of the survival probability decreases with 
increasing a. 

All the exponents of the generic transition do not de- 
pend on p or a, implying that the critical behavior is 
universal. Moreover, their values are in excellent agree- 
ment with the known values of the PC (or DP2) uni- 



D. Transition at a = 

After discussing the generic transition, we now turn to 
the line of transitions at pf < p < p* and a = 0. To 
investigate these transitions more closely, we performed 
both spreading and density decay simulations at fixed p 
and several a- values approaching a — (as indicated by 
the solid (red) arrows in the phase diagram, Fig. [I}. 

Let us start by discussing the density decay simula- 
tions. Figure [7] shows the stationary density p s t of active 
sites as a function of a for several values of the heal- 
ing rate p. Interestingly, the stationary density depends 
linearly on a for all healing rates p^, p < p < p,*, in 
seeming agreement with mean-field theory. This means 
p st — B^<t u with uj — 1 and B^ being a /i-dependent 
constant. We also analyzed, how the prefactor B^ of the 
mean-field-like behavior depends on the distance from 
the simple contact process critical point. As inset a) of 
Fig.[7]shows, B^ diverges as (p — p^, p )~ K with k = 2.3(1). 

At the critical healing rate p cp of the simple con- 
tact process, the stationary density displays a weaker 
cr-dependence. A fit to a power-law p st ~ a Wc ' p gives 
an exponent value of uj cp = 0.108(2). In contrast, at the 
endpoint at healing rate p* , the corresponding exponent 
u>* = 1.4(1) is larger than 1. 

These results of the density decay simulations must be 
contrasted with those of the spreading simulations. Fig- 
ure[S]shows the time dependence of the survival probabil- 
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FIG. 7: (Color online) Density decay simulations. Main 
panel: stationary density p 3t as a function of the boundary 
rate a for various healing rates p. For /x° p < p < p*, the solid 
lines are fits of the low-er behavior to p„t = B^a. At the sim- 
ple contact process critical point, p — p^ p = 0.30324, and at 
the endpoint, p = p* = 0.552, we fit to power-laws p st ~ cr" 
which gives exponents of a; C p = 0.108(2) and uj* = 1.4(1). 
The data are averages over 50 to 200 runs with system sizes 
L = 2000 to 5000. Inset a: prefactor of the linear a de- 
pendence as a function of p — p cp . A fit to a power law gives 
B M ~ (p-p c ,?y K with k = 2.32(10). Inset b: prefactor B M as 
a function of p* — p. A fit to a power law gives B M ~ (p* —p) K 
with k* = 0.91. 



ity P s for p — 0.4 and several a. At early times, all curves 
follow the <7 = data due to the small values of the rate 
of the boundary activation process (j4}. (Note that the 
(7 = curve does not reproduce the survival probabil- 
ity of the simple contact process. This is because in our 
generalized contact process, a sample is surviving as long 
as not every site is in state Ii even if there are no active 
sites.) In the long-time limit, the P s curves approach 
nonzero constants, as expected in an active phase. How- 
ever, in contrast to the stationary density p st (Fig. [7]), the 
stationary value of P s does not go to zero with vanish- 
ing boundary a. Instead, it approaches a tr-independent 
constant. We performed similar sets of simulations at 
other values of p in the range p^ p < p < p* , with anal- 
ogous results. We therefore conclude that - somewhat 
surprisingly - the survival probability and the station- 
ary density of active sites display qualitatively different 
behavior at the a = phase transition. 

We now show that the properties of these quantities 
can be understood within a simple domain wall theory. 
The relevant long-time degrees of freedom at p > pf 
and a <C 1 are the domain walls between Ii and I2 do- 
mains. These domains are formed during the early time 
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FIG. 8: (Color online) Spreading simulations: Survival prob- 
ability P 3 as a function of time t at p = 0.4 for various values 
of the boundary rate a. The data are averages over 100000 
runs. Inset: Low-cr limit of the stationary P 3 as a function of 
p. The dashed line is a fit to P s ~ (p* — p)^ with p* = 0.552 
and j3 = 0.87(5) in agreement with the PC universality class 
(see, e.g., Refs. HH). 



evolution when the system follows the simple contact pro- 
cess dynamical rules ([T]) to ([3]). At late times, the do- 
main walls can hop, they can branch (one wall branching 
into three), and they can annihilate (two walls vanish 
if the meet on the same bond between two sites). This 
means, the domain wall dynamics follows the branching- 
annihilating random walk with two offspring (BARW2). 

In our case, the BARW2 dynamics is controlled by two 
rates, the domain wall hopping rate T and the branching 
rate fl (annihilation occurs with certainty if two walls 
meet). These two rates depend on the underlying gen- 
eralized contact process dynamics. In the limit a <^ 1 
they are both linear in the boundary rate, T = aFr{p), 
Q = (tFq,{p) because a single boundary activation event 
is sufficient to start a domain wall hop or branching (Fr 
and F(i are nontrivial functions of p). Because both rates 
are linear in er, their ratio is a- independent, thus the 
steady state of the domain walls does not depend on a in 
the limit a < 1. This explains why the survival proba- 
bility P s of the generalized contact process saturates at a 
nonzero, <7- independent value in Fig. [8l It also explains 
the cr-dependence of the stationary density p s t of active 
sites in the following way: For a <C 1 and p > ac- 
tive sites are created mostly at the domain walls at rate 
<7. Consequently, their stationary density is proportional 
to both cr and the stationary domain wall density pa-wi 
i.e., p s t ~ crpdw, in agreement with Fig. [7J (The linear 
CT-dependence of p s t is thus not due to the validity of 
mean- field theory.) 

These results imply that the phase transition line at 
a = between p^P and p* is not a true critical line be- 
cause there is no (nontrivial) diverging length scale. It 
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only appears critical because the stationary density of ac- 
tive sites vanishes with a. Note that this is also reflected 
in the fact that the system is not behaving like a criti- 
cal system right on the phase transition line a — (no 
power-law time dependencies, for instance). Instead, the 
physics of this transition line is controlled by the BARW2 
dynamics of the domain walls with a finite correlation 
length for all p? < p < p* . 

E. Scaling at the contact process critical point 

(Mc p ,0) 

Even though the generalized contact process is not crit- 
ical at a = and p > p?, its behavior close to the critical 
point of the simple contact process can be understood in 
terms of a phenomenological scaling theory. 

Let us assume that the stationary density of active sites 
close to (p?,0) fulfills the homogeneity relation 

p st (Afx,a) = b^^pstiAfib- 1 ^,^-^) (10) 

where Ap = p — p° p and b denotes an arbitrary scale 
factor. /3 cp and v~. are the usual order parameter and 
correlation length exponents and y cp denotes the scale 
dimension of a at this critical point. Setting b — a 1 / ycp 
then gives rise to the scaling form 

p st (Ap,a) = aP»p/(.<v«p)x (Ana-UK*™^ (11) 

where X is a scaling function. At criticality, A/i = 0, this 
leads to p st (0,cr) - cr<WK>cp) ( using X (0) = const). 
Thus, uj cp — Pcp/i^cpVcp)- For a — > at nonzero A/i, 
we need the large-argument limit of the scaling function 
X. On the active side of the critical point, A/i < 0, 
the scaling function must behave as X(x) ~ \x\^ cp to 
reproduce the correct critical behavior of the density, 
p st ~ \p- IX? _ 

More interesting is the behavior on the inactive side of 
the critical point, i.e., for A/i > and a — > 0. Here, we 
assume the scaling function to behave as X(x) ~ x~ K . 
In this limit, we thus obtain p st ~ (Ap)~ K a u (just as 
observed in Fig. [7]) with oj = ((3 cp + K)j{v^ v y cp ). As 
a result of our scaling theory, the exponents ui,ui C p and 
k are not independent, they need to fulfill the relation 
w cp (/3 cp + k) = ficpUj. Our numerical values, uj = 1, 
uj cp = 0.108(2) and k = 2.32(10) fulfill this relation in 
very good approximation, indicating that they represent 
asymptotic exponents and validating the homogeneity re- 
lation (TO}. Using p cp = 0.2765 and v^ p = 1.097 [H|, the 
resulting value for the scale dimension y cp of a at the 
simple contact process critical point is y cp = 2.34(4). 

F. The endpoint (/AO) 

Finally, we turn to the point (//*, a) = (0.552, 0) where 
the generic phase transition line terminates on the /i axis. 



At first glance, one might suspect this point to be a mul- 
ticritical point because it is located at the intersection of 
two phase transition lines. However, we argued in Sec. 
IIVDI (based on the domain wall theory) that the transi- 
tion line at a = and p cp < p < p* is not critical. This 
implies that the endpoint (ji* , 0) is not multicritical but 
a simple critical point in the same universality class (viz., 
the PC class) as the generic transition at p > p* . In fact, 
the endpoint can be understood as the critical point of 
the BARW2 domain wall dynamics in the limit a — > 0. 

To test this hypothesis, we first study the survival 
probability and density of active sites as p* is approached 
along the p axis. The inset of Fig. [8] shows the stationary 
survival probability (more precisely, its saturation value 
for a — > 0) as a function of p. The data can be well fitted 
by a power-law P s — (p*-p) 13 with f3 = 0.87(5). The cor- 
responding information on the stationary density of ac- 
tive sites can be obtained from inset b) of Fig. [7] It shows 
the prefactor £? M of the linear tr-dependence p s t = B^a 
as a function of p* — p. Sufficiently close to p*, their 
relation can be fitted by a power law Bp ~ (p* — p) K 
with k* — 0.91. Thus both /3 and n* agree with the order 
parameter exponent of the PC universality class within 
their error bars. This confirms the validity of the domain 
wall theory of Sec. II V Dl at p* . 

The discussion of the cr-dependence of P s and p s t right 
at p* is somewhat more complicated because it is deter- 
mined by the subleading er-dependencies of the domain- 
wall rates T and ft. Moreover, because the dynamics 
is extremely slow at p w p* and u < 1, our numeri- 
cal results close to the endpoint are less accurate then 
our other results. According to the domain wall theory 
of Sec. IIV Dl the stationary survival probability should 
fulfill the homogeneity relation 

P s (Ap,<r) = bV ,/± P s (Apb- 1 / l ' ± ,ab-v*) (12) 

where A/i = p — p* while j3 and v 1 - are the order pa- 
rameter and correlation length exponents of the BARW2 
transition (PC universality class). The only unknown ex- 
ponent is y* . The same homogeneity relation should hold 
for the domain wall density, but not the density of active 
sites. 

Setting the scale factor to b = a x l y gives the scaling 
form 

P s (Ap,a) = a^^y'^YiApa-^^y'^) . (13) 

Right at the endpoint, Ap = 0, this gives P s - a ^/^ ±v '\ 
To test this power-law relation and to determine y* , we 
performed spreading simulations at p = p* and several 
a between 0.03 and 1. The low-tr behavior (not shown) 
can indeed be fitted by a power law in a with an expo- 
nent /3/(u y*) = 0.5(1). Using the well-known values 
/3 = 0.92 and = 1.83 of the PC universality class, we 
conclude y* = 1.0(2). Within the domain wall theory, 
Pdw ~ -Ps and the stationary density of active sites is 
Pst - o~ Pdw - ^" with oj* = 1 + /3/(^V) = 1-5(1). 
This agrees well with the numerical estimate of 1.4(1) 
obtained from the density decay simulations in Fig. [7J 
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The scaling form (fT3)) can also be used to determine 
the shape of the phase boundary at [i > [i*. The phase 
boundary corresponds to a singularity of the scaling func- 
tion Y at some nonzero value of its argument. Thus, the 
phase boundary follows the power law a ~ (/i — [i*Y y . 
At fit of the data in Fig. pleads to v^y* = 1.8(2) which 
implies y* = 1.0(1) in agreement with the above estimate 
from the spreading simulation data. 

To investigate the time dependence of P 8 close to the 
endpoint, the homogeneity relation (1121) can be general- 
ized to include a time argument. On the right hand side, 
it appears in the scaling combination (t/to)b z with to the 
basic microscopic time scale. It is important to realize 
that this microscopic scale diverges as cr _1 with a — > 
(independent of any criticality at fj,*). Thus, the right 
scaling combination is actually tab 7 - . We used the result- 
ing scaling theory to discuss the power-law decay of P s on 
the phase boundary shown in Fig. [5^- The scaling the- 
ory predicts P s ~ <r~^t~ s with ( = 5 as the endpoint is 
approached. This agrees with our numerical data (shown 
in the inset of Fig. [5^) which give £ s» 0.284 

In summary, all our simulation data support the no- 
tion that the endpoint (/z* , 0) is a not a true multicritical 
point but a simple critical point in the same universal- 
ity class (PC) as the entire generic phase boundary at 
ft> fj,*. The behavior of some observables makes it ap- 
pear multicritical, though, because the microscopic time 
scale of the domain wall dynamics diverges with a —> 0. 



V. CONCLUSIONS 

In summary, we have studied the phase transitions of 
the generalized contact process with two absorbing states 
in one space dimension by means of large-scale Monte- 
Carlo simulations. We have found that this model has 
two different nonequilibrium phase transitions, (i) the 
generic transition occurring for sufficiently hight values 
fj, > fx* of the healing rate and nonzero values of the 
boundary activation rate a, and (ii) a transition at ex- 
actly a — for [if < fi < fi* . 

The generic transition is in the parity-conserving (PC) 
universality class (which coincides with the DP2 class in 
one dimension) everywhere on the /i > u* phase bound- 
ary, in agreement with earlier work [l3l. Il7j. In contrast, 
the cr = transition turned out to be not critical. The 
density of active sites rather goes to zero with the vanish- 
ing boundary activation rate a while the survival prob- 
ability remains finite for a — > 0. Its behavior is con- 
trolled by the BARW2 dynamics of the domain walls be- 
tween different inactive domains (which is not critical for 
[if < /i < (J,*). It is interesting to note that the be- 
havior of our model at a = differs qualitatively from 
the o — y limit of the finite-a behavior in the entire 
parameter region [if < fj, < fx*. 

As a result, the crossover between directed percolation 
(DP) critical behavior at a = and parity conserving 
(PC) critical behavior for a > does not take the naively 



expected simple scaling form. In particular, the generic 
(a > 0) phase boundary does not continuously connect 
to the critical point of the a = theory (the simple con- 
tact process critical point). Instead, it terminates at a 
separate endpoint (/z*,0) on the /z-axis. While this point 
shares some characteristics with a multicritical point, it 
is actually just a simple critical point in the same univer- 
sality class (PC) as the entire generic phase boundary. 

We emphasize that the crossover between the DP and 
PC universality classes as a function of a in our model is 
very different from that investigated by Odor and Meny- 
hard (20| . These authors started from the PC universality 
class and introduced perturbations that destroy the sym- 
metry between the absorbing states or destroy the parity 
conservation in branching and annihilating random walk 
models. They found more conventional behavior that can 
be described in terms of crossover scaling. In contrast, 
the transition rates (p} to (fj} of our model do not break 
the symmetry between the two inactive states anywhere 
in parameter space. 

Crossovers between various universality classes of ab- 
sorbing state transitions have also been investigated by 
Park and Park (21 '2'A . They found a discontinuous 
jump in the phase boundary similar to ours along the 
so-called excitatory route from infinitely many absorbing 
states to a single absorbing state [2l[. Moreover, there is 
some similarity between our mechanism and the so-called 
channel route [22j from the PC universality class to the 
DP class which involves an infinite number of absorbing 
states characterized by an auxiliary density. In our case, 
at a = (but not at any finite a), any configuration 
consisting of Ii and I2 only can be considered absorbing 
because active sites cannot be created. The density of 
I1-I2 domain walls then plays the role of the auxiliary 
density; it vanishes at the endpoint (/i*,0). However, 
our crossover occurs in the opposite direction than that 
of Ref. [22|: The small parameter a takes the system 
from the DP universality class to the PC class. Note 
that an unexpected survival of active sites has also been 
observed in a version of the nonequilibrium kinetic Ising 
model with strong disorder. Here, the disorder can com- 
pletely segment the system, and in odd-parity segments 
residual particles cannot decay [24| ■ 

The generalized contact process as defined in eqs. (P) 
to (j4|) is characterized by three independent rates (one 
rate can be set to one by rescaling the time unit). In the 
bulk of our paper, we have focused on the case ft = [i for 
which our system reduces to the usual contact process in 
the limit of a — > 0. In order to study how general our 
results are, we have performed a few simulation runs for 
ft =/= [i focusing on the fate of the endpoint that sepa- 
rates the generic transition from the a = transition. 
The results of these runs are summarized in Fig. |9] which 
shows the phase diagram projected on the ft — ft plane. 
The figure shows that the line of endpoints of the generic 
phase boundary remains distinct from the simple contact 
process (tr = 0) critical line in the entire ft— ft plane. The 
two lines only merge at the point ft — 0, fi — 1 where the 
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FIG. 9: (Color online) Projection of the phase diagram of the 
generalized contact process on the jl — fi plane. The individual 
symbols show the locations of the phase boundaries as deter- 
mined from our simulations: solid blue circles - transition for 
(j = (simple contact process), solid red triangles - generic 
transition for a = 1, open squares - approximate location of 
the endpoint of the generic transition (a — ► 0) estimated from 
the transition at a = 0.01. The lines are guides to the eye 
only. Points A and B are the simple contact process critical 
point and the endpoint investigated in the main part of the 
paper. 



system behaves as compact directed percolation |13j ]. 

Our study was started because simulations at /j fj% p 
and it < 1 [18[ seemed to suggest that the generalized 
contact process with two absorbing states is always active 
for any nonzero a. The detailed work reported in this 
paper shows that this is not the case; a true inactive 
phase appears, but only at significantly higher /i > fi*. 
Motivated by this result, we also carefully reinvestigated 
the generalized contact process with n — 3 absorbing 
states which has been reported to be always active (for 
any nonzero a) in the literature [l3l Il7j |. However, in 
contrast to the two-absorbing-states case, we could not 
find any inactive phase in this system. 

Let us close by posing the question of whether a sim- 
ilar splitting between the n = 1 critical point and the 
ri = 2 phase transition line also occurs in other micro- 
scopic models with several absorbing states. Answering 
this questions remains a task for the future. 
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